Light scattering in inhomogeneous Tomonaga-Luttinger liquids 
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We derive the dynamical structure factor for an inhomogeneous Tomonaga-Luttinger hquid as can 
be formed in a confined strongly interacting one-dimensional gas. In view of current experimental 
progress in the field, we provide a simple analytic expression for the light-scattering cross section, 
requiring only the knowledge of the density dependence of the ground-state energy as they can be 
extracted e.g. from exact or Quantum Monte Carlo techniques, and a Thomas-Fermi description. 
We apply the result to the case of one-dimensional quantum bosonic gases with dipolar interaction 
in a harmonic trap, using an energy functional deduced from Quantum Monte Carlo computations. 
We find an universal scaling behavior peculiar of the Tomonaga-Luttinger liquid, a signature that 
can be eventually probed by Bragg spectroscopy in experimental realizations of such systems. 



I. INTRODUCTION 



It is well known theoretically that systems of re- 
duced dimensionality, especially in one dimension, 
present simultaneously enhanced quantum fluctuations 
and stronger interaction effects that can lead to exotic 
ground states [l], [2[. ^From the experimental point of 
view, there are many prototypical one-dimensional _^ys- 
tems, that range from organic 0, 0] or inorganic 
conductors, and antiferromagnetic (AF) spin chain 
or ladder 0, [l^] materials, to nanoscale systems such 
as of quantum wires [III i carbon nanotubes [l3l - [l6| 
or self organized An atomic wires on Ge(OOl) semicon- 
ductor surfaces [l3l. More recently, advances in atom 
trapping technology has permitted the realization of both 
fermionic and bosonic one-dimensional systems with un- 
precedented control [lB-[2l|. The low-energy physics of 
such one-dimcnsional systems is well described by the 
Tomonaga-Luttinger Liquid (TLL) theory [i[i,[ll[ll. 
In a single component TLL, there is a single gapless 
branch of excitations with linear dispersion, and the in- 
terplay between interactions and quantum fluctuations in 
the ground state leads to power-law decay of correlations 
with interaction dependent exponents. Remarkably, the 
low energy theory is fully characterized by two parame- 
ters: the velocity u of the linearly dispersion excitations, 
and the dimensionless exponent K controlling the decay 
of all correlations, the corresponding exponents being ra- 
tional functions of K . In physical systems several promi- 
nent features of TLL have been observed after measuring 
the spectral function 0, |^ , the structure factor [H or the 
conductivity [l3| . and more recently the first quantita- 
tive check of TLL physics has appeared for the spin- 1/2 
ladder material bis(pipcridinium) tetrabromocuprate(II) 



(C5Hi2N)2CuBr4 (abbreviated BPCB in the following), 
in an applied magnetic field [l^. However, despite this 
recent achievement, in many of the physical systems men- 
tioned above, little control can be exerted on the values 
of u and K and thus the Luttinger exponent K is taken 
as an adjustable parameter [^, Il3| . This fact prompts 
for the search of more than one signature of Tomonaga- 
Luttinger liquid physics for a single system. 

In the case of systems with strong confinement (e.g. 
confined quantum gases), excitation properties can be 
most easily accessed by light spectroscopy techniques, as 
proposed in the early days of atomic Bose-Einstein Con- 
densation [13, d^. For example, the spectral function 
has recently been measured in trapped Fermi gases by ra- 
diofrequency spectroscopy [l^ and the dynamical struc- 
ture factor has been studied successfully by optical Bragg 
spectroscopy in free and trapped Bose-Einstein conden- 
sates [27l - [3]| as well as trapped Fermi gases [12] . Bragg 
spectroscopy can be based on energy transfer to the sys- 
tem at fixed momenta (ssl - fsst or can permit the study of 
the full momentum composition of excitations by a coher- 
ent momentum transfer mapping [s^. For these reasons, 
Bragg spectroscopy can be especially useful to investi- 
gate the properties of the many phases realizable in these 
systems such as Mott insulator, Tonks-Girardeau gas 
or supersolid phases as recently proposed [37l - [45| . The 
most recent experimental progress in producing long- 
lived ground-state polar molecules in a three-dimensional 
(3D) optical lattice, and possibly also in ID arrays of pan- 
cakes and 2D arrays of tubes [i^ as well as condensates 
of dipolar atoms j47l |48| , opens up wide perspectives in 
the comprehension of controlled quantum systems with 
tunable short and long range interactions under progres- 
sively reduced dimensionality. 
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As we have more extensively reviewed in Ref. [491 ] . 
among many possible realizations, quantum dipolar gases 
in ID confinement are quite peculiar TLL systems. Here 
in fact, one single parameter drives the crossover from 
weak to strong interaction regimes, where however the 
weakest regime is a Tonks- Girardeau state, the strongest 
being a Density Dipolar Wave state characterized by 
quasi-ordering. Based on the above motivations, we de- 
rive an analytic expression for light-scattering intensity 
in the case of a weakly inhomogeneous TLL. This expres- 
sion is valid within a Thomas-Fermi description, where 
the system can be considered locally homogeneous. The 
expression requires the knowledge of the density depen- 
dence of the ground-state energy of the homogeneous sys- 
tem, as can be obtained by e.g. approximate calculations, 
exact Bethc-Ansatz technique or Quantum Monte Carlo 
(QMC) simulations. The paper is organized as follows. 
After reviewing in Sec.|ll]the calculation of the dynamic 
structure factor and the inelastic light-scattering cross 
section of homogeneous Tomonaga-Luttinger Liquids, we 
derive in Sec. IIIII the general expression for the inhomo- 
geneous system within the Thomas-Fermi approach, in 
terms of the eigenvalues and eigenfunctions of the hydro- 
dynamic TLL. We then specialize in Sec. IIVI to the case 
of one-dimensional quantum bosonic gases with dipolar 
interaction in an harmonic trap, using our previous QMC 
findings [i^. Here the results are explicitly discussed in 
the various regime while the single parameter built-up 
from density and interaction strength is tuned. 



Earlier experimental studies [28| have shown that con- 
densate properties of atomic cold gases could be studied 
by means of Bragg scattering yielding large energy res- 
olution and sensitivity. The system is illuminated by 
two lasers beams of momenta ki and ^2 and frequencies 
LJi , UJ2 of difference uj that creates a periodic field whose 
intensity is proportional to cos [(fci — k2) ■ r — ujt]. The 
external potential couples to the density n{q) of the sys- 
tem where q — ki — fc2- After using the golden rule, the 
response of the system to this perturbation is the dynam- 
ical structure factor [5T| . Light scattering experiments 
then directly measure S{q,uj). 

This quantity is then a benchmark against the theo- 
retical descriptions of the systems. For an homogeneous 
Tomonaga-Luttinger liquid occuring in interacting one- 
dimensional system the dynamic structure factor can be 
readily obtained [l^. In the following, we briefly sketch 
the derivation. For a system of interacting spinless parti- 
cles, either bosons or fermions, the low-energy physics is 
that of a Tomonaga-Luttinger liquid whose Hamiltonian 



, dx 

""-It. 



A 



(3) 



with u the velocity of the excitations and K the 
Tomonaga-Luttinger exponent. The density operator 
n(x) is expressed in terms of bosonic operators (j): 



n(x) = no dxip + A,n coa{2m((f>(x) ~ Trnox)), (4) 



II. LIGHT-SCATTERING CROSS SECTION IN 
HOMOGENEOUS TOMONAGA-LUTTINGER 
LIQUIDS 



The dynamic structure factor S{q,u;) is central in the 
description of interacting many-body systems. S{q, uj) 
is related to the Fourier transform of the imaginary 
density-density correlation function with the fluctuation- 
dissipation theorem. It is therefore accessible by means 
of inelastic scattering, where density fluctuations are in- 
duced in the system and their subsequent relaxation is 
measured revealing the system characteristics. While in- 
elastic neutron scattering has been the tool to probe the 
condensate nature of superfluid helium and the roton 
spectrum js^], inelastic light scattering has been pro- 
posed and widely used in dilute quantum degenerate 
gases. Within linear response theory the scattering cross 
section a of light at frequency uj and angle il incident on 
a Bose atomic sample is: 

oc —{nsiuj) + l)Imx(g,w) = S{q,uj) , (1) 

where nsi^j) is the Bose distribution function, n = N/V 
and x((7, lo) is the Fourier transform of the density-density 
correlation function 



with m an integer and the equilibrium density. 

If the wavelength of the incoming light is much larger 
than the average interparticle distance, we can neglect 
the contribution of the oscillatory terms in Eq. Q . Using 
translational invariance, the expression for the density- 
density response function becomes: 

Xix - x\ t) = t),dx4>{x\ 0)]). (5) 

Knowing that the time-ordered correlation function 
(r,[(/)(x,T) -(/)(0,0)]2) = KFi{x,t) with Fi[x,t) = 
log[(.T^ -I- (u|t| -I- a)^)/a^]/2, the imaginary part of the 
response function ^ can be obtained [l^l as 



lmx{q.oj)^^uK[5{io + u\q\)-5{io-u\q\)] , (6) 
giving the scattered intensity at zero temperature: 



cc S{q,uj) = sign(aj)Imx((j', cj) 



K\q\ 



[Sioj + u\q\) + Sioj - u\q\)]{7) 



X{r,t) ^ ~i0it){[n{r,t),n{O,O)]). 



(2) 



Expression ([7]) embodies the symmetry with respect to 
inversion of the velocity u as required by Galilean in- 
variance, and evidences the dependence of the light- 
scattering signal from the ratio q/uj. 
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III. LIGHT-SCATTERING CROSS SECTION IN 
INHOMOGENEOUS TOMONAGA-LUTTINGER 
LIQUIDS 



A. Hydrodynamic approach 

The presence of an external potential V{x) confining 
the cold atomic cloud induces density inhomogeneity, 
and the external light perturbation probing the density- 
density correlation function introduces time-dependent 
processes. The treatment of the problem is easier under 
conditions of weak inhomogeneity and slow processes as 
they can be met in experiments, where external poten- 
tials vary on length and time scales longer than the char- 
acteristic system quantities, and local equilibrium hydro- 
dynamic behavior sets in. Under these conditions, the 
gas can be still described by a hydrodynamic Tomonaga- 
Luttingcr Liquid Hamiltonian [i^, [52H56j 



_ dx 

J^T LL — I — 

Zir 



u{x) 
K{x) 



{d^cj,{x)f 



(8) 

Here, the boundary conditions imposed arc (j){~R) = 
and = —ttN, with N the number of particles in 

the system. The parameters u{x) and K(x) now de- 
pend on position. In analogy with the homogeneous case, 
where u and K are related by the expressions u/K = 
(hTr)^^d^/dn and by Galilean invariance uK = Trhn/m, 
one sets: 



u{x)K{x) 

u{x) 
K{x) 



TT—no[x) 
m 

1 ( d[i{n) 



flTT \ d 



(9) 
(10) 



n=no(x) 



Once an estimate of the equilibrium density tiq [x) and of 
the chemical potential ^{n) are known, this phenomeno- 
logical approach allows the determination of u{x) and 
K{x). 

The response function ([S]) in the case of the Hamilto- 
nian ([8]) can be calculated using the decomposition: 



(\dx'^^^ 

J— it u{x') 



N 



J-R""^ u{x') 



(11) 



Here, [a„,aJJ = dn,m and the first term comes from the 
addition of N particles in the system. The functions (p„ 
satisfy the eigenvalue equation: 



- Lolifn = u{x)K{x)dx 



u(x) 
K{x) 



(12) 



with boundary conditions (/5„(±i?) = 0, and the normal- 
ization 



V3„(x)(p,„(a;) 

dx . . . = On-. 



(13) 



u(x)K(x) 

The influence of the trapping potential enters eq. (jl2p via 
the equations for u(x) and K{x) ([9]). The density-density 
response function thus can be expressed as: 

X(x, x', t) ^ 9{t) ^ ^^^sinM, (14) 
^-^ nojn dx dx' 

Taking the Fourier transforms with respect to x and 
x' and the Laplace transform with respect to t, we find: 



27raj„ 



1 



1 



(15) 



where Im(z) > 0. Finally, taking the limit z — > cj + iO+ 
we obtain: 



lmxiq,UJ + iO+) ^ |- ^ |^„((3')|^[(5(w - w„) 



2uj 

n 



(16) 



Eq. (|16p maintains the structure of its homogeneous 
counterpart ([6]). 

The density-density response function can be deter- 
mined whenever the density dependence of the ground 
state energy per unit length e{n) or of the chemical po- 
tential fJ.(n) = {§f^)\n=n(x) is known. An especially sim- 
ple situation is realized when e{n) oc rp~^^. That type of 
dependence of energy on density corresponds to several 
limiting cases of ID TLL systems. For example, in the 
Lieb-Liniger gas [HI, [III there are two well understood 
limits. At low density or strong repulsion, the gas be- 
haves as a hard-core boson gas [59| with 7 = 1, while at 
high density or weak repulsion, the Bogoliubov approxi- 
mation applies and gives and energy density proportional 
to n^, so that 7 = 0. The study of the crossover between 
these two limits requires the Betlie-Ansatz computation 
of the ground state energy density [H^l- A similar situa- 
tion occurs in the case of dipolar gases. For low densities, 
the energy per unit length e(n) has the 7 = 1 behav- 
ior typical of hard core bosons, while for high density it 
has the 7 = 2 behavior of a crystal of classical dipoles, 
and a Dipolar-Density- Wave manifests [i^. As density 
increases, the system crosses over from the low density 
hardcore boson gas to the high density Dipolar-Density- 
Wave. 

In the model with e{n) = grO^'^ and in the case of har- 
monic trapping potential V{x) = mQQX^ /2, the eigen- 
values uJn of (|1H) can be found exactly, and the func- 
tions ip„ are expressible in terms of Gegenbauer polyno- 
mials [59, roOl as: 



fnix) 



2 \ a + 1/2 



i?2 



\R 



1) 



(17) 
(18) 
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Here, uq and Kq are the Tomonaga-Luttinger parameters 
corresponding to the density at the trap center, 



^" = \/ 7:T{n + 2a + 2) ^ ^ r(l + a) , 



(19) 



and a = (7+1)^^ — 1/2. In particular, in the case of hard- 
core Bose gas when 7 = 1, a = and the Gegenbauer 
polynomials reduce to Chebyshev polynomials [6l|. In 
order to calculate the scattered light intensity, we need 
the Fourier transform of the (finS. Using Eq. (7.321) of 
Ref. ^ we obtain: 



\f>n{q)\^ = 2uoKoR{n + a + 1) 



r(n + 2a + 2) 

r(n + i) 



{qR) 



2a+2 



(20) 



where the Jm ai'c the Besscl functions of the first kind. 
Thus: 

uqKq v-^ , , r(n + 2q; + 2) 

lmx(g,w + iO+) = -^^(n + a + 1) ^ ' 



Roj 



Jn+a+l{qR) 



[qR]'- 



[5{U} - lUn) + 5{lO + LOn)] 



r(n + 1) 

(21) 



Eq. ((2T|) shows the main features of the scattered light 
intensity. This is a set of discrete peaks, whose weight 
is a function of gi?, and whose spacing reduces with in- 
creasing the trap size R ^ 00. 



B. Approach via Density- Functional Theory with 
Local Density Approximation 

In the present section we derive an approximate ex- 
pression for the dynamical structure factor of an inho- 
mogeneous ID TLL, reverting to the Density Functional 
Theory (DFT) accompanied by a Local Density Approx- 
imation (LDA) . We sketch in the following the main con- 
cepts and derivation. Through the Hohenberg and Kohn 
theorem, DFT establishes that the ground state energy 
of a system subjected to an external potential V{x) is a 



hmctional Eg[n{x)] 



E[n{x)] + n{x)Vix)dx of the 
density n{x), where E[n(x)] embodies the kinetic and 
exchange-correlation parts. The equilibrium density pro- 
file is determined by the variational condition 



5Eg[n{x)] 
Sn{x) 



(22) 



stating that equilibrium corresponds to a minimum of the 
energy against changes in the particle density, while the 
total number of particles is fixed through the (density- 
dependent) chemical potential ^. Eq. (|22p reminds the 
Thomas-Fermi equilibrium condition in non-interacting 
systems, and in fact the Density Functional sets a one- 
to-one correspondence between the ground state energies 



of an interacting system and of its non-interacting ana- 
logue. Whenever an analytic expression of n{n) is avail- 
able, inversion of the equation of state ([^^ allows the 
determination of the equilibrium density ng (x) ■ 

While Eq. ([2^ is exact, the actual determination 
of the functional E[n{x)] needs approximations. Under 
the conditions of shallow confinement, we can safely use 
the Local Density Approximation. Here, the functional 
E[n{x)] is replaced by 



E^^^[n{x)] 



[n{x)]n{x)dx 



(23) 



where e'"""(ri,) is the energy per particle of the homoge- 
neous system with density n. 

Differentiating £;^^^[n(a;)] = Eg[n{x)] + J dx{V{x) - 
X)n{x) with respect to n{x), A being a Lagrange multi- 
plier fixing the total number of particles, one obtains the 
condition for the local chemical potential 



^i[n{x)] ^ V{R) - V{x), 



(24) 



where the local chemical potential is defined by the func- 
tional derivative: 



SE 
6n(x) 



dij 



^hotn 



in)) 



dn 



(25) 



n—no{x) 



If an analytic expression of /x(n) is given, Ea. (|25p would 
allow to find n{x) by inverting the relation Eq. ([M)) . The 
energy e^°'^{n) can be obtained after perturbation the- 
ory, or by exact calculations such as Bethe-Ansatz, or 
else by computational Quantum Monte Carlo methods. 

We now turn to the problem of determining the dynam- 
ical structure factor of the inhomogeneous system. To 
this aim, we follow the reasoning in psl [63| and imag- 
ine to slice it into small segments of length Ax, where 
the density no(x) can be considered uniform, and thus 
sum together all the contributions ([7]) of the different 
segments. The dynamical structure factor of the inho- 
mogeneous system would then be approximated by: 



M-.^-.^ni.x)) 



(26) 



S\ioia{qi'-^in) is given by Eq. ([7]), where now the 
Tomonaga-Luttinger parameters u = u(n) and K = 
K{n) depend on density. 

With the help of ((71), we obtain: 



S{q,uj) 



\q\ 



dxK {nQ(x)) 



[5{lj - u{n(3{x))\q\) + 5{lo -t- w(no(a;))|g|) 



(27) 



Introducing x*{Lj/\q\)^ such that lo = u[n[x*))\q\ we can 
rewrite: 



S{q,uj) 



K{no{x*)) 



1 du 1 


1 d7io 1 


1 dn In— no (a;*) 


1 dx \ . 



(28) 
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Since the compressibility is a positive quantity, the chem- 
ical potential is an increasing function of the density. 
Moreover for a trapping potential that is an increasing 
function of position, from Eq. 1^^) the density is seen 
to decrease with position. Thus, when the velocity is an 
increasing function of density, the solution x* turns out 
to be unique 

The quantity can be obtained by differentiating 



(29) 



the relation (1241) with respect to x , i.e.: 



/ n=no(x) 

We can therefore write 



dno 
dx 



dV 
dx 



K'inoix*)) 



d^ 



n— no (x*) 



1 du 1 


dV 


1 dn \ n—nQ{x*) 


\ dx \ 



(30) 



We now use the relation u^in) — —4-4 obtained from 
^ and rewrite ([50)) as: 



S{q,uj) 



Trh 



no{x*) 



R 



— I 

dx \x=x'' 



1 -|_ 



(31) 



with the notations e'{x) = de/dn, e"{n) = d^ejd'n?, and 
e"'(n) = d'^e/dit^. 

Formula pip represents the main result of this paper. 
It gives an analytical expression for the light scattering 
cross-section of an inhomogeneous TLL once the ground 
state energy as a function of the density is known, e.g. 
by an exact analytical (Bethe-Ansatz) or via numerical 
simulations (QMC). Remarkably, Eq. ([51]) predicts that 
S{q,uj) is only a function of w/lqj. In fact, this is the spe- 
cific signature of Tomonaga-Luttinger Liquid behavior in 
shallow trapped ID Bose systems, as it can be measured 
by Bragg spectroscopy. 

In order to illustrate the relevant features and make the 
connection with Eq. (j2ip obtained via the hydrodynamic 
approach of Sec. lIII Al we now treat the case of harmonic 
trapping. In this case dV/dx = milpX, and using Eq. ([M)) 
we have mfll\x*\ = \/2mQl{e' {no{0)) — e'(p*)), where 
we have set p* — no{x*) and u{p*) ~ ^/\q\- Eq- ([5T|) 
thus simplifies into: 



S{q,uj) 



nhp* 



R^2mnl{e'inom - e'ip*)) 



'(P* 



{32) 



'(p. 



We now check the consistency of the result ([5T|) with 
2T|) . by explicitly calculating (|32|) for the model e{Ti) oc 



n7 + 2 



. Eq. 



then reads: 



S{q,uj) 



(- 



(7+l)mr!2i?2 V25(7 + 2) 



n„(0)=10 



no{0)=l 




25 

q (a.u.) 50 



FIG. 1: (Color online). TLL model with e(nro) oc rf and 
7 = 2 in a harmonic trap. S{q, to) in arbitrary units in the 
(tj,g) plane and different densities at the trap center. 



that at the trap center. We first notice that the dynami- 
cal structure factor in p3p makes explicit the characteris- 
tic already embodied in the structure of Eq. ([32)) . namely 
that S{q, uj) depends on wavevector and frequency solely 
through their ratio w/g. Second, the formula (j33|) with 
7=1 agrees with the result of Ref. [g^l , in the limiting 
ijj 3> 9^/2 case. Finally, in App. |^we show by inspection 
that the LDA approximation ([55)) is fully recovered from 
expression ([^. 



Fig.|T]displays the 3D plot of 5(g, w) resulting from the 
use of ([55]) in the (w, q) plane, while varying the densities 
at the trap center. S{q, uj) is a set of discrete peaks whose 
position varies linearly with lj j q and such linear behavior 
is independent on the interaction strength. 



1 



V"o9 ) 



where we have defined uq = u(no(0)) as the velocity of 
excitations in a uniform system having a density equal to 



Before proceeding to apply Eq. ([321) to a dipolar ID 
Bose gas, we step on commenting the found correspon- 
dence between hydrodynamic and DFT-LDA approaches 
on a more general footing. It is well known for nor- 
mal Fermi systems [g^] with extension to Bose super- 
fluids [65| . that the treatment of dynamical processes in 
interacting inhomogeneous systems do require the devel- 
opment the Current-Density Functional Theory, where 
invariance conditions render the energy to be a func- 
tional of the current besides density. It was demon- 
strated that the analogue of LDA leads in this case to 
Navier-Stokes equations (Landau-Khalatnikov two-fluid 
equations for superfluids), where viscosities, densities and 
currents (normal and superfluid) have a microscopic ex- 
pression in terms of Kubo relations and low-frequency 
response functions as they can be calculated in the ho- 
mogeneous system at the local densities and currents. 
Such a general view is reflected by the present result. 
In the Tomonaga-Luttingcr-Liquid free harmonic Hamil- 
tonian, where the interactions are effectively embodied 
in u and K , the Navier-Stokes equations become indeed 
the simple hydrodynamic relations of Sec. IIII Al On the 
other hand, in the DFT and LDA approach of Sec. IIII Bl 
the treatment explicitly uses the two mappings: from 
interacting to non-interacting system (DFT) and from 
inhomogeneous to homogeneous (LDA). 
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IV. ID BOSE GASES COUPLED VIA DIPOLAR 
INTERACTIONS 



In this Section we specialize to the case of a ID dipolar 
gas in a harmonic trapping potential. We first recall the 
main results known for the homogeneous system, and 
then apply Eq. ([5^ to determine the scattered light 
intensity. The system is characterized by the strength 
of the interactions Cddi resulting from cither magnetic 
Cdd = A*oA*d O'" electric Cdd = d? /^o dipolcs, where ^d and 
d are the magnetic and electric dipole moments and /xq 
and eo are the vacuum permittivities. An effective Bohr 
radius can be defined from Cdd as rg = MCdd/ C^tt^"^) 
and the Hamiltonian in effective Rydberg units Ry* = 
n^/{2Mrl) is 



H = {nrof 



i<3 



, (34) 



where lenghts are expressed in 1/n units. The physics of 
the model is entirely specified by the dimensionless cou- 
pling parameter nrg, so that in the high-density limit the 
system becomes strongly correlated and a quasi-ordered 
state occurs, where the potential energy dominates. 

The ground-state energy e(n) of this model was deter- 
mined by means of Reptation QMC method in Ref. [i^ . 
In the low nro — > limit it reproduces the Tonks- 
Girardeau (TG) state energy per particle of a free spin- 
less Fermi gas, whose energy per particle is eTcin) = 
7r^(nro)^/3 Ry* ■ In the large nrQ — > oo limit of high- 
density dipolcs, it reproduces the Dipolar Density Wave 
(DDW) state where eoDwin) = C(3)(nro)^ Ry* and 
C(3) = 1.20205. The QMC thermodynamic energy per 
particle in Rydberg units can be represented as an ana- 
lytical function of nr^: 



ep{nro) 



C{3){nro)'^ + a{nro) 



1 -I- nro 



bjnro y + c{nro)(^+a^^ 



3 l + d{nro)9 



The fitting coefficients, yielding a reduced — 5, are: 
a = 3.1(1), b = 3.2(2), c = 4.3(4), d = 1.7(1), e = 
3.503(4), / = 3.05(5), and g = 0.34(4). 

The Bragg intensity is thus easily obtained by Eq. ((32)) 
once the value of p* is determined. In Figl^] we report 
the scaling behavior of S{q,Lu) vs. uj/{uoq) for differ- 
ent densities at the trap center uq. Larger rigro indi- 
cate stronger coupling interactions, crossing over from 
TG to DDW states. The linear behavior in the low 
Lo/q regime is striking, the slope continuously increas- 
ing with decreasing u(0) and thus uq. In the TG limit, 
the tail of S{q,uj) is insensitive to changes of the den- 
sity at the center of the trap, and in fact the curves with 
'^o'^o = 0-01 and 0.1 do coincide. The comparison with 
the TG gas (7 = 1 and norg = 0.01) and the DDW case 




0.4 0.5 0.6 

o)/[u(0) q] 

FIG. 2: (Color online). ID dipolar Bose gas confined in an 
harmonic trap, with e(?i) as determined by QMC simulations. 
S{q,ij) (arbitrary units) vs. uj/{uoq) for different densities no 
at the trap center. The values of no'^o running from 0.01 to 
1000 are indicated in the legend. 



(7 = 2 and rtoro = 10'^) is better seen in FiglJl where 
5'((j, Lo) 1/ e' (tt-o)/ no is plotted as a function of Lu/{uQq). 
One can notice that a crossover takes place in the in- 
termediate densities regime. Viewed in the log-log scale, 
the plot evidences how a measure of the S{q,uj) tail to- 
wards small Lu/q, would provide a way to determine the 
interaction regime. A peculiarity of the TLL behavior 
is the power-law trend when u/iu^q = 1) is approached. 
A detailed study of the power-law non-analyticity for a 
trapped Bose gas can be found in [45| . 



100 



3 
w 




0.0001 



0.0001 



0.001 0.01 
(o/[u(0)q] 



FIG. 3: The same as Fig. [21 but in a log-log scale. The com- 
parison with the Tonks-Girardeau limit gas (7 = 1) and the 
dense dipole limit corresponding to a Dipolar Density Wave 
(7 = 2) is shown in evident manner. 
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V. CONCLUSIONS 

Wc have derived the dynamical structure factor for 
an inhomogeneous Tomonaga-Luttinger liquid as it can 
occur in a confined strongly interacting one-dimensional 
gas. In view of current experimental progress in the field, 
we have provided an easy-to-use and simple analytical 
expression for the light-scattering cross section, Eq. ([31]l , 
valid within a Local Density Approximation. 

The analytical expression ([31]) predicts that S{q,uj) is 
only a function of uj/\q\ and is the central result of this 
work. In fact, this is the specific signature of Tomonaga- 
Luttinger Liquid behavior in shallow trapped ID Bose 
systems, along with a power-law behavior when oj /{u^q) 
is approached, as it can be measured by Bragg spec- 
troscopy. 

Expression (j3ip is validated by the independent deriva- 
tion (|2ip by means of a hydrodynamic approach, which is 
reported in detail in App.lA] The connection between the 
two approaches is a second result of this work, and is a 
consequence of the more general Current-Density Func- 
tional Theory [g^, applied to the conditions of the 
present work. 

We thus remark that expression (j3ip can be in principle 
applied to the many ID systems cited in the introductory 
material, once the trapping potential is known together 
with the ground state energy as a function of the density, 
e.g. by means of perturbative, exact, or computational 
methods applied to the homogeneous system. Extension 
of the present method to include additional local pertur- 
bations coupling to the density, could be used to investi- 
gate the propagation of local density fluctuations. 

Finally, we have applied our findings to the case of 
one-dimensional quantum bosonic gases with dipolar in- 
teractions, using the harmonic profile typical of experi- 
ments in this field, accompanied by our previous QMC 
data for the energy per particle. We find an universal 
scaling behavior peculiar of the Tomonaga-Luttinger liq- 
uid [401, a signature that can be eventually probed by 
Bragg spectroscopy in ongoing experimental realizations 
of such systems [4q . 



Appendix A: Justification of the LDA formula 

In order to justify the approximate formulas, it is more 
convenient to work with the integrated intensity: 



i{q,n) 



S{q, uj)dLO 



(Al) 



since the delta functions in the sum ([2T|) contribute as 
step functions in I{q, il) giving more regular expressions. 
Using the approximation ([55]). we expect: 



/(g,r!) 



Trhuonolql 2a + I 



n 

uoq 



2(a + l) 



X2F1 



1 / n 

a + l,-;a + 2; 

2 \UQq 



(A2) 



where 2F1 is the Gauss hypergeometric function. Using 
the expression ([2T|). we obtain instead the exact expres- 
sion: 



I{q,ujN) oc > ^ '-^ A3 



n<N 



where N is such that wat = f2. In order to check the 
consistency between (|A2[) and (|A3p we can work on the 
sums in (|A3[) . We expect that the sum is dominate d by 



the terms having n ^ 1. Using Ec^. (9.3.2) in Rcf. [61 
we expect that for qR < n, Jn+a+iilR) is an exponen- 
tially small quantity with n. For qR > n however, Eq. 
(9.3.3) in Ref. [6l| suggests that: 



(A4) 



7r(ri -1- a + 1) tan/3 
X cos^((n + a + l)(tan/3 - /3) - 7r/4) , 

where cos /S = qR/ {n + a + 1). Elementary trigonometry 
gives the approximation: 

{n + a + l) 



TTy^iqRy - (n + a + l)2 
X 1 1 + sin 



(A5) 



— (n + a + 1) arccos 



2v/(gi?)2-(n + a + l)2 

n + Q' + 1 



qR 



Dropping the term oscillating with n in ([A5|) . we use as 
approximation: 

0{qR -n-a-1) 



r 



1(9^) 



We can also approximate: 



(A6) 



r{n + 2a + 2) 

r(n+ 1) 
n + a + I 



n 



2a+l 



yj{n + l)n + 2q; + 1 
So that the sum in Eq. (|A3[) can be approximated by 



(A7) 
(A8) 



E 

n<min(Af,g_R) 



t2q+1 



yJ{qRY -{n + a + lf 



(A9) 



Finally, by approximating the sum ([A9p by an integral, 
we find: 



E 



,2a+l 



^{qRf -{n + a + lf 



n<Tam(N ,qR) 

nmm(7V,gii) ^^^2^+1 



(AlO) 
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Using Eq. (6.6.1) in [6i[, we have (for qR < N): 

/O 



^0 



1 



(gi?)2"+iB(j,/,^,,^,).(a + 1,1/2), 



(All) 



where Bj.{a,b) are the incomplete Beta functions. With 
Eq. (6.6.8) of [6l| we can check that Eq. ()A2p agrees with 



the obtained approximate expression (|A11[) . 
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